The Discovery of 7-ray Emission from the Blazar RGB J0710+591 
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The high-frequency-peaked BL Lacertae object RGB J0710+591 was observed in 
the very-high-energy (VHE; E > 100 GeV) wave band by the VERITAS array of 
atmospheric Cherenkov telescopes. The observations, taken between December 2008 
and March 2009 and totahng 22.1 hours, yield the discovery of VHE gamma rays 
from the source. RGB J0710-I-591 is detected at a statistical significance of 5.5 stan- 
dard deviations (5.5cr) above the background, corresponding to an integral fiux of 
3.9 ± 0.8 X 10^12 g-i ^2,% of the Crab Nebula's flux) above 300 GeV. The ob- 

served spectrum can be flt by a power law from 0.31 to 4.6 TeV with a photon spectral 
index of 2.69 ± 0.26stat i 0.20sys. These data are complemented by contemporaneous 
multiwavelength data from the Fermi Large Area Telescope, the Swift X-ray Telescope, 
the Swift Ultra- Violet and Optical Telescope and the Michigan-Dartmouth-MIT ob- 
servatory. Modeling the broad-band spectral energy distribution with an equilibrium 
synchrotron self-Compton model yields a good statistical fit to the data. The addition 
of an external-Compton component to the model does not improve the fit nor bring 
the system closer to equipartition. The combined Fermi and VERITAS data constrain 
the properties of the high-energy emission component of the source over four orders 
of magnitude and give measurements of the rising and falling sections of the spectral 
energy distribution. 

Subject headings: gamma rays: galaxies — BL Lacertae objects: individual (RGB J0710-I-591, 
EXO 0706.1+5913, VER J0710+591) 



1. Introduction 



RGBJ0710+591 was originally dis covered by HEAP A -1 (IWood et al.lll984l l: it is a BL Lac 
object located at a redshift oi z = 0.125 (|Giommi et al.lll99lh . Hubble Space Telescope observations 
of this blazar show a fully re solved diskless gala xy that can be modeled using an elliptical profile 
with a nuclear point source (jScarpa et al.ll200d ). In addition, RGBJ0710+591 is found to have 
an optical 'compact companion' located at a distance of 1.4" (3.5 kpc) from the core of the active 
galactic nucleus (AGN). The nature of this companion is unknown, but it could be a foreground 



star o r an unresolved faint gala xy. It is, however, to o luminous for a globular cluster (jFalomo et al 
2000) and is not seen in radio (iGiroletti et al.ll2004l ). 



Blazars are characterized by spectral energy distributions (SEDs) consisting of two peaks, the 
low energy peak arising from synchrotron emission and the high energy peak fror n either leptoni c 
or hadronic interactions (for a review of the emission mechanisms in blazars, see iBottcherl 120071 ) . 
Depending on the location of the low energy peak, BL Lacs are classified as being high-frequency- 
peaked (HBL, fpeafe ~ 10^^^^^ Hz), intermediate-frequency-peaked (IBL, Vpeak ~ 10^^^^^ Hz), or 
low-frequency-peaked (LBL, Vpeak ~ 10^^~^^ Hz). The location of the low energy SED peak of 
RGBJ0710+591 (identified as EXO 0706.1+5913 in the following reference) is higher than ~ lO^*^ 



-7- 



Hz (jNieppola et al.ll2006l ). clearly identifying it as an HBL. 



RGB J0710+591 is a well known HBL featured in many catalogs, but it was not detected in 
the high-energy (HE; E > 100 MeV) band by the EGRET instrument on the Compton Gamma- 
Ray Observatory. In the very-high-energy (VHE; E > 100 GeV) band, the Whipple 10-m imaging 
atmospheric Cherenkov telescope (lACT) reported several flux upper limits for observations taken 
in different epochs: (0.91, 1.69, 3.58, and 4.29) xlO~^^ cm~^ s~^ (99.9% confidence level) above 
(350, 350, 500 and 400 GeV) respectively , corresponding to (8.7%, 16.1%, 52.4% and 45.9% ) of 
the Crab Nebula's flux (jHoran et al.ll2004l ). The HEGRA lACT system reported a flux upper limit 



above 1.08 Te V of 0.91 x 10 cm ^ s ^ (99% confidence level) corresponding to 6% of the Crab 



Nebula's flux (jAharonian et al.ll2004l ) 



There are several compelling reasons to study RGB J0710+591 in the MeV to TeV energy 
range. Due to the hardness of its X-ray spectrum, the high energy peak is ex pected to be located 
in an ideal location for HE and VHE observations. Th ere are both lept o nic (iGhisellini &: Madau 
19961 : bottcher fc ChianeJ bood ) and hadronic models (jAharonianl boool : iMiicke et al] bood ) that 
explain the photon-production mechanism of the high energy peak. The continuous measurement 
of the SED from MeV to TeV energies is vital to constrain certain aspects of the emission models 
and to perhaps exclude some of them entirely. For example, some leptonic models appear harder 
at HEs than hadronic models. At VHE energies, photons intera ct via pair-production w ith the 
infrared component of the extragalactic background light (EBL) (iGould &: Schrederl 119671 ). This 
causes a decrease in the observed flux and a softening of the observed spectrum. Thus, extragalactic 
objects at moderate redshifts with reasonably hard VHE spectra are ideal candidates to indirectly 
study the history of t he EBL and its SED, particularly in the optical to far infrared (see, e.g., 
Aharonian et al.l 120071 ). 



VHE 7-ray emissi on from RGBJ0 710+591 was discovered by VERITAS during the 2008- 
2009 observing season (lOng et al.ll2009l ). Following this, the VERITAS collaboration initiated a 
multiwavelength (MWL) observation campaign that included the participation of the Fermi Large 
Area Telescope (LAT) in HE 7-rays, the Swift X-ray Telescope (XRT) in X-rays and the Swift Ultra- 
violet and Optical Telescope (UVOT) and the Michigan-Dartmouth-MIT (MDM) observatory in 
optical. The resulting MWL SED spans thirteen decades in energy. 



VERITAS Observations 



VERITA S is an array of four 12-m lACTs loca ted in southern Arizona at the F. L. Whipple 
Observatory (Iweekes et ahlbood : iHolder et ahlbood ). VERITAS has an energy resolution of ~15% 
between 100 GeV and 30 TeV, an angular resolution of ~0.1° (68% containment) per event, and 
a field of view (FoV) of 3.5°. VERITAS can detect a source with a fiux of 1% of the Crab in < 50 
hours, while a 5% Crab Nebula flux source is detected in ~2.5 hours. These speciflcations are valid 
for the data taken during this campaign; the VERITAS array was subsequently upgraded, resulting 
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in a sig: 



;nificant improvement in sensitivity (jPerkins et al.ll2009l ). 



RGBJ0710+591 was observed with VERITAS from December 2008 to March 2009 dm'ing 
astronomical darkness and partial moonlight for a total time of 34.9 hours. Approximately 30% 
of these observations are rejected due to poor weather or hardware-related problems with the 
majority rejected due to weather, resulting in 22.1 hours of exposure (livetime corrected). The 
average zenith angle of the observa tions is 30°. The VERITAS data are calibrated and cleaned 
as described in I Daniel et al.l (|2008l ). After calibration, several noise-reducing selection criteria 
(cuts) are made, including rejecting those events where only the two closest-spaced telescopes 
participated in the trigger . The VER ITAS standard analysis consists of stereo parametrization 
using a moment analysis jHillas 1985) and following this , scaled parameters are used for event 
selection ( Aharonian et al. 19971 : Krawczynski et al. 2006 ). This selection is based on a set of 
cuts that are a priori optimized using data taken on the Crab Nebula where the 7-ray excess 
has been scaled by 0.01 to account for weaker sources like RGB J0710+591. VERITAS observed 
RGB J0710+591 in "wobble" mode, where the center of the FoV is offset by 0.5° fr om the blazar, 
to allow for si multaneous source and background measurements (jBerge et al.l 120071 ) The reflected 



region model (jAharonian et al.ll200ll ) is used to calculate the background events from several off- 
source regions. A total of 576 on-source events are observed from the direction of RGB J0710+591 
and 3890 are seen in the background regions. Depending on the pointing direction either eight 
or nine background regions are used resulting in a final a of 0.115 (a is the ratio of the areas 
of the on and off-source regions). This results in an excess of 129 events. This excess has a 
statistical significance of 5.5 standard deviations {a) and corresponds to an integral flux of 3.9 it 
0.8 X 10-12 s-i (3% 

of the Crab Nebula's flux) above 300 GeV. There is no evidence for 
variability during these observations (a fit of the lightcurve to a constant flux yields a reduced 
chi-squared value of 0.21). The excess is point-hke and is located at 07^" 10"" 26.4" ± 2.4^^^^, 
59° 09' 00" ± 36stat ( J2000, see Figure EJ the systematic uncertainty in the pointing is 90") which 
is 49" away from the location of RGB J0710+591 (07^" 10"^ 30.0 7^59° 08' 20.5") measured by 
the Very Large Array at 5 GHz ([Laurent-Muehleisen et al.l Il999l ) . The VERITAS source name 
is VER J0710+591. The distributions for the on-source and background regions are shown in 
Figure [1] {0 is the distance between the source location and the location of an event) . The photon 
energy spectrum can be fit by a simple power- law {dN/dE = I{E/Eq)-'^, Eq = I TeV) with photon 
index Fvhe = 2.69 it 0.26stat i 0.20sys (see Tabled] for the spectral points). The fit yields a per 



degree of freedom of /dof 



Jsys 

1.3/3. 



Fermi Observations 



The Fermi LAT is a pair-conversion gamma-r ay detector seiisitive to photons in the energy 
range from below 20 MeV to more than 300 GeV (jAtwood et ahlbood ). The LAT data for this 
analysis were taken between 4 August 2008 and 4 August 2009. The data are analyzed using the 
standard likelihood tools distributed with the Science Tools v9rl5p3 package available from the 
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Fig. 1. — The distribution of 9^ for on-sourcc VHE events (point (shaded region) from observations 
of RGB J0710+591. 9 is the distance from the source location to the event location. The dashed 
line represents the location of the cut on 6^ applied to the data. The 9'^ cut retains events which 
lie to the left of the line. 



Table 1: The VERITAS spectral bins. 'E;o^' and 'Ehigh axe the low and high energy bounds, 'Flux' 
is the differential flux, 'Excess' is the number of excess events and 'Sig.' is the significance of the 
bin. 
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Right Ascension [deg] 

Fig. 2. — Smoothed VHE significance map of RGB J0710+591 from the VERITAS observations. 
The VLA location of the blazar is shown as a star. The solid circle is centered at the best fit 
location of the excess emission in the VERITAS data; the radius of the circle is the sum of the 
statistical and systematic uncertainties. The dashed circle is the best fit plus a 95% containment 
radius of 2.7' to the Fermi LAT data. Both locations are consistent with the VLA position of 
RGB J0710+591. The VERITAS point spread function is shown in the top right corner. 
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Fermi Science Support Center 0. Only events having the highest probabiUty of being photons, 
those in the "diffuse" class, are used. To limit contamination by Earth albedo gamma rays which 
are produced by cosmic rays interacting with the upper atmosphere, only events with zenith angles 
< 105° are selected. Events with energy between 100 MeV and 300 GeV and within a circular 
region of 10° radius centered on the V LA coordinates of RGB J0710 +591 are analyzed with an 



unbinned maximum likelihood method ( Cash 



197S 



Mattox et al .11 19961 ). The background emission 
is modeled using a Galactic diffuse emission model and an isotropic component □ The fluxes are 
determined using the post-launch instrument response functions P6_V3_DIFFUSE. 

A point source is detected with a significance of more than 8 standard deviations. The best-fit 
position (07^ 10°^ 37% 59° 09' 32") has a 95% containment radius of 2.7' and is statistically consistent 
with the VLA coordinates of RGB J0710-I-591 (see FigureE]). The highest energy photon associated 
with the source has an energy of 75 GeV. The time-averaged Fermi spectrum is fit by a power-law 
for which / = 1.43 ± 0.35stat ± O.lOsys x 10^^^ cm'^s"^ MeV-\rHE = 1-46 ± 0.17stat ± 0.05sys, 
and Eq = 8.8 GeV is the energy at which the correlation bet ween the fit value s of I and The is 
minimized. The spectral index reported in the IFGL catalog (jAbdo et al.ll2010l ) for this source is 
1.28 ± 0.21stat- The discrepancy between that result and the result reported here can be explained 
by the different energy range considered in the likelihood model of the source: 100 MeV - 100 GeV 
for the catalog analysis, 100 MeV - 300 GeV for the present analysis. The lack of photons above 75 
GeV carries information that is important for the fit results. No evidence for variability is detected. 



4. Swift Observations 



A total of seven Swift (jGehrels et al.l 12004] ^ XRT and UVOT observations of RGB J0710-h591 
were performed between 20 February 2009 and 2 March 2009 (observation identifications 00031356001 
through 00031356007). The total XRT observation time is 15.2 ks while the total UVOT obser- 
vation time is 14.8 ks. The Swift XRT data are analyzed with HEAsoft 6.5 using the most recent 
calibration files. The first three XRT observations in photon-counting mode have source count 
rates of ~0.7-1.5 counts s~^, resulting in significant photon pile-up in the core of the point-spread 
function, which occurs for rates >0.5 counts s~^. The photon pile-up is avoided by using an annular 
region with inner radius ranging from 2 to 7 pixels and an outer radius of 30 pixels (47.2") to extract 
the source counts. Background counts are extracted from a 40 pixel radius circle in a source- free 
region. For the last four XRT observations in windowed-timing mode, source counts are extracted 
from a rectangular region of 40 by 20 pixels high, with background counts extracted from a nearby 
source-free rectangular region of equivalent size. Evidence for marginal flux variability between 
the Swift XRT observations is indicated by a probability of 0.004 for a fit of a constant fiux. The 
time-averaged XRT data can be fit with an absorbed power-law model where the hydrogen column 



""^ http://fcrmi.gsfc. nasa.gov/ssc/ 

^http://fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels. html 
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,20 2 ( Kalberla et al.l l2005l ) . resulting in a combined spectrum 



density Nh is fixed at 4.4 x 10 
with photon index of 1.86 it 0.01, with a x^/dof = 371 / 303. This hard spectrum is consistent 
with a synchrotron peak above 10 keV, as can be inferred from the radio to X-ray s pectral energy 



distr ibution (SED) generated from archival data, including an EXOSAT detection ([Giommi et al, 
I995I ). 



Swift UVOT (jRoming et al.ll2005l ) data are presented in the six bands of v, b, u, uvml, uvm2 
and uvw2. The UVOTSOURCE tool is used to extract counts, correct fo r coincidence losses, apply 
background subtraction, and calculate the source flux (jPoole et al.ll2008l ). The standard 5" radius 
source aperture is use d, with a 20^^ backgro und region. The source fluxes are extinction corrected 
using the dust maps of lSchlegel et al.l (jl998l ). Since the stacked UVOT images are not deep enough 
to sample the host galaxy profile, no attempt is made to correct for the contribution of the host 
galaxy and this can be seen as an increase in fiux in the optical portion of the SED (Figure [3]) 
starting with the u-band. 



MDM Observations 



The observations of RGB J0710+591 taken from February 19 through 24, 2009 with the MDM 
Observatory 1.3m telescop^ in R and B band are bias-subtract ed and fiat-fielded using the Im- 
age Reduction and Analysis Facility (IRAF) (|Todvl 119931 . Il986l ) . Instrumental magnitudes are 
extracted using the IRAF DAOPHOT package and calibrated using four comparison stars in the 
field, with compari son star magnitudes f rom the Naval Observatory Merged Astrometric Dataset 
(NOMAD) catalog (Izacharias et al.ll2004l ). and extinction coefficients Ab = 0.166, = 0.106 from 
the NASA/IPAC Extragalactic Database (NEDfl- 

The host galaxy of RGB J0710+591 makes a significant contribution to the optical fiuxes 
measured by MDM. The host galaxy subtractions for this BL Lac were calculated previously in 
Hvvonen et al.l (j2007l ) and these subtractions are used to correct for the host galaxy contribution 
to the fiux measurements. The resulting extinction- and host-galaxy-corrected optical fiux in the 
R-band (4.68 x 10^"^ Hz) is 3.51 ± 0.07 x 10^^ Jy Hz and in the B-band (6.85 x 10^"^ Hz) is 4.23 ± 
0.41 x 10^^ JyHz. Even with the host galaxy corrections, the MDM data are still higher than the 
UV measurements made by the Swift UVOT which are dominated by the nucleus. Thus, the MDM 
and Swift optical data in the u, b and v bands still have contributions due to the host galaxy and 
should be considered upper limits. No substantial optical variability is seen over the course of the 
five day observation period. 



http://www.astro.lsa.umich.edu/obs/mdm/ 
^http://ned www.ipac.caltech.edu/ 
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Fig. 3. — Broad-band SED of RGB J0710-I-591. The closed squares are R- and B-band data from 
the MDM 1.3m telescope, the open triangles are from the Swift UVOT, closed circles are from 
the Swift XRT, the shaded area is the Fermi LAT measurement and the open circles are the 
measured VERITAS data . The broad-band spectral model is the equilibrium SSC version from 
Bottcher &: Chiang (120021 ) with a cooling break energy occurring around 1.3 x 10^*^ Hz. Note that 
the optical measurements (MDM and Swift u, b and v bands) are affected by emission from the host 
galaxy and a re thus not used in t he SE D fit. The SED model curve includes EBL absorption using 
the model of iFranceschini et al.l (|2008l ) and the dashed line is the de-absorbed VHE spectrum. 
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6. Modeling and Discussion 



The observed broad-band SED of RGB J0710+591 is shown in Figure [3l The data, excluding 
the MDM and Swift u, b and v bands which are affected by contributions from the host, are fit with 
a purely lept onic synchrotron se l f-Com pton (SSC) mod el. Specifically, an equilibrium version of 
the model of Bottcher &i Chiang ( 2002 ). as described in lAcciari et alJ ( 2009 ) is used. The best-fit 
parameters are listed in Table [21 The VHE data in the plot are the measured va lues and the model 
i nclud es the effects of intergalactic 7 — 7 absorption using the EBL models of iFranceschini et al. 
(|2008l ) and the reported redshift of z = 0.125. De-absorbing the VHE data using the same models 
results in a de-absorbed spectral index of Fvhe ~ 1.85. 

Many of the model parameters are not constrained by the observations. In particular, the 
choice of the bulk Lorentz factor (Fbuik) and observing angle (^o6s) are somewhat arbitrary and 
similarly good fits can be obtained using other values. Evidence from portions of the radio structure 
is supportive of a small observing angle which corroborates the s mall value obtained h ere. However, 
the bulk Lorentz factor derived from these same data is lower (jGiroletti et al.ll2006l ). Regardless, 
the synchrotron portion of the SED always turns over far below the X-ray regime. Therefore, in 
order to obtain the observed hard X-ray spectrum, a very hard injection index (q = 1.5) has to be 



assum ed. This poses challenges to standard (Ist-order) Fermi acceleration models (jStecker et al 



20071 ). and it might indicate second-order Fermi acceleration or other modes of acceleration (e.g., 
shear at boundary layers). The fit requires a rather low magnetic field of B = 0.036G , far below 
equipartition. This cannot be remedied even with the choice of a higher Doppler factor. 

Additionally, an external-Compton (EC) component was added to the model to see if the 
resulting fit would improve, or if the system would be closer to equipartition. However, in order 
to produce an EC peak at sufficiently high frequencies, a very low magnetic field is needed. Thus, 
the inclusion of an EC component does not lead to an overall improvement of the fit statistically, 
or with respect to equipartition. 

With a photon index Fvhe = 2.69, RGB J0710+591 is one of the hardest VHE blazars de- 
tected to date. Furthermore, out of the 19 blazars d etected by both V HE i nstruments and Fermi , 
RGBJ0710+591 has the hardest HE photon index (|Abdo et ahlbood ) and iNieppola et al.l kood ) 
found that RGB J0710+591 has one of the highest low energy SED peaks out of a sample of over 
300 BL Lacs and place it in a list of 22 ultra-high-frequency-peaked candidates. Nevertheless, the 
observed SED is well fit by a purely leptonic SSC model without the need to include an external 
component. 

Observations of th e similarly hard (Fvhe = 2.50) and distant {z = 0.1396) blazar lES 0229+200 
(jAharonian et al.l 120071 ) indicated that the EBL density in the mid-IR band was close to the lower 
limits measured by Spitzer. These constraints on the EBL are based on the assumption that 
the intrinsic spectrum of blazars must have a spectral index of 1.5 or larger. Thus, if blazars like 
lES 0229+200 and RGB J0710+591 are unusual objects this assumption might not be valid because 
their intrinsic spectral indices could be smaller than 1.5. However, the modeling of RGB J0710+591 
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and other similar blazars indicates that hard-spectrum VHE blazars might not be exotic (i.e. 
having an intrinsic spectral index less than 1.5) and are part of the general blazar population. 
Therefore, using VHE observations of hard spectrum blazars to constrain EBL models as done in 



Aharonian et al.l (|2007l ) continues to be a valid undertaking. 
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Table 2: The equilibrium model parameters for the purely leptonic SSC model that was used to fit 
the SED of RGB J0170+ 591. Details of the model can be found in the tex t. 

Parameter Value 



7mm 6 X 10^ 

Imax 2 X 10^ 

eT Injection spectral index 1.5 

Escape time parameter r/esc = 100^ 

Magnetic field at 0.036 G 

Bulk Lorcntz factor Tbuik = 30 

Blob radius 2 x 10^^ cm 

Qohs 1.91*^ degrees 

Redshift z = 0.125 

Le (jet) 4.49 X lO''^ erg/s 

Lb (jet) 1.75 x 10^^ gj.g/s 

Ls/Le 3.90 x 10-2 



^equivalent to the superluminal angle 



